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ABSTRACT 

The recent localization of some short-hard gamma ray bursts (GRBs) in galaxies with low star 
formation rates has lent support to the suggestion that these events result from compact object binary 
mergers. We discuss how new simulations in general relativity are helping to identify the central engine 
of short-hard GRBs. Motivated by our latest relativistic black hole-neutron star merger calculations, 
we discuss a scenario in which these events may trigger short-hard GRBs, and compare this model 
to competing relativistic models involving binary neutron star mergers and the delayed collapse of 
hypermassive neutron stars. Distinguishing features of these models may help guide future GRB and 
gravitational wave observations to identify the nature of the sources. 

Subject headings: black hole physics — gamma rays: bursts — gravitational waves — stars: neutron 



> 

m 
o 

o 

Oh: 

6 



1. INTRODUCTION AND OBSERVATIONS 

Gamma-ray bursts are short duration phenomena 
characterized by time- varying, high-energy, non-thermal 
electromagnetic emission. Based on their duration 
and energy spectra they are typically categorized 
into two categories: "short /hard" and "long/soft" 
ijKouveliotou et alJl99.'^ . Many optical and X-ray coun- 
terparts of long bursts have been seen, some in co- 
incidence with Type Ib/c supernovae, repre senting ex- 
tremely en ergetic collapses of massive stars (|Paczvnskil 
Il998t iHiorth .2003.') . By contrast, it has been much 
more difficult to identify counterparts for short GRBs 
(SGRBs), except for the known soft gamma repeaters, 
which may be observed as SGRBs (Palmer et al. 200^, 
but cannot explain more th an a smal l fractio n of the ob- 
served SGRB sample ( Lazzati et al.l 12005: ..Nakar et al.l 
1200513). 

Recently, the Swift and HETE-2 satellites have local- 
ized for the first time the X-ray aft erglow following short- 
period GRBs: G RBs 05 0509b (Gehrels et al.' '200l, 
050724 ( Barthclmv et al.l [2005: Berger et al. .2000). 
0508 13 lIBerged I2005ir and 0512 21 JSoderberg et all 
mi^ by Swift; GRB 050709 l|Hiorth et abl I200,5t 
Fox et al. 2005) by HETE-2. Details about the physi- 



cal parameters observed and inferred from these bursts 
can be found in larger (2006). In all cases, the inferred 
host had a rather low star-formation rate. This find- 
ing disfavors the collapse of a massive star as a progen- 
itor, since those systems have very short lifetimes. In- 
stead, it favors the identification of a compact object 
bi nary me rger as the progenitor, as originally suggested 
by iPaczvn ski (1986); a significant fraction of both neu- 
tron star-neutron star (NSNS) and black hole-neutron 
star (BHNS) binaries will take longer than 1 Gyr be - 
tween formation and merger l|Belczvnski et al.l I200I . 
Other scenarios, such as the accretion-induced collapse 
of a NS with a WD companion have been proposed 
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ijDermer fc Atov;iiill2006D . but have not been studied nu- 
merically in as much detail. 

The observed characteristics of SGRBs can be exam- 
ined in light of population synthesis calculations for com- 
pact object binaries, but it is difficult to disentangle 
BHNS vs. NSNS merger scenarios given the current ob- 
servational and theoretical uncertainties. The inferred 
rate of observed SGRBs and the merger rate of their 
putative sources are consistent, but the latter are un- 
certain by at least two orders of magnitude: TZgrb = 
(r2/47r)7?.,„ = (0^/2)7?.„i, where 9 is the beaming angle. 
The most recent GRB rate estim ates yield TtnRR y 1 — 
3 Myr-^ per Milky Way Galaxy ijGuetta fc Piranll2006l: 
INakar et al J 1200531) . over an order of magnitude large r 
than previous results l|Guetta fc P iranl2005iH Andcl2004() . 
Breaks in the observed x-ray spec tra imply ~ 0.2 — 0.3, 
and thus VI/At: = 0.02-0.05 ("Sod erberg et al. 2006^. The 
predicted compact object merger rates per Milky Way 
Galaxy faU in the range 7^,„ - 1 - 100 Myr'^ fo r NSNS 
and n,n ~ 0.5 - 10 M yr-^ fo r BHNS iVoss fc TaurisI 
12003: Belc zvnski et 31.1 12002: Ka logera et al.ll2004f ). The 
rates are uncertain by at least an order of magnitude and 
the beaming angle by a factor of a few. An extremely 
high SGRB rate may indicate that at least some SGRBs 
are the products of NSNS mergers, as these are much 
more common tha n BHNS mergers in population syn- 
thesis calculations ()Voss fc Tauri£l2003llBelczvnski et al] 

Hoi. 

The projected distance from a localized SGRB to the 
center of its host galaxy is insensitive to whether the 
merging binary is an NSNS or BHNS system. Recent re- 
sults show the respective projected distance likelihood 
functions are nearly the same for large galaxies, and 
have only minor differences in smaller ones, which would 
rule out p ositive identi fications based on location alone 
l|Belczvns ki et al.ll200(J) . 

The range of fluxes seen from short GRBs is rather 
large, indicating that the burst energy must depend 
sensitively on at least one parameter of the progenitor 
system. The observed fluxes of short GRBs typically 
fall in the range 5 x lO^'' — 10~^ erg cm~^ s~^. an d 
fluences lO'^ '^ - lO^^ '^ erg cm'^ (|Balazs et al■ll2003^ . 
Assuming a characteristic distance of 1 Gpc, we see 
that the luminosities L and total energy released E for 
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the burst satisfy (47r/0)L = 10 



50.5 



10 erg s and 



(47r/rj)£' = lO'^^-^ - 10^°-^ erg, respectively. 

The first locahzed short burst, GRB 050509b, has 
an extremely low measured isotropic energy in 7-rays, 
= 3 X {n / 4:Tt) erg, compared to previously ob- 
served short GRBs. Assuming that the GRB is at the 
measured redshift of z = 0.225, the luminosity of the 
burst is actually similar to the other localized bursts, 
~ 10^° erg s~^, but it lasted for a m uch shorter t i me. In 
general, as we see from Table 1 of lJanka et all (p005), 
there is significantly smaller variation in the peak lumi- 
nosity than in the isotropic energy output. 

We note that the low energy for GRB 050509b ar- 
gues for an extremely low density of baryons surrounding 
the GRB. Typically, it is assumed that the Lorentz fac- 
tor of the jet r will be no bigger than the ratio rj of 
the energy in the jet to the rest energy of the baryons 
through which it mu st travel, i.e. T < rj = E^/MbC^ 
l|Shemi fc Piranll990() . For GRB 050509b, the least ener- 
getic of the observed bursts, assuming the Lorentz factor 
r > 100 (Occhslin & Janka 2005) implies that at most 
lO~*(Sl/47r)M0 of material can surround the progenitor. 
T hese results are co nfirmed by the numerical calculations 
of lAlov et al.1 1)2005(1 ■ who find that the density of baryons 
surrounding the burst must fall rapidly with increasing 
radius. For sufficiently large baryonic loading, they find 
the observational counterpart to a merger is not a GRB. 

2. PROGENITOR MODELS; THEORETICAL AND 
OBSERVATIONAL CONSTRAINTS 

Standard models for merger-induced short-duration 
GRBs involve hot, massive accreti on disks (M > 
O.OIM©) around spinning BHs (see IPiranI 120051 for a 
thorough review). As the disk is accreted, one of two 
possible mechanisms is responsible for the creation of a 
gamma-ray jet. One suggestion is that the hot material 
emits neutrino-antineutrino pairs that annihilate above 
and below the disk to produce a relativistic jet co ntain- 
ing e~ — e"*" pairs and photons (see, e.g. |Paczvnsk]| ll986t 
IGoodman et al.iri987|) . Calculations of accreting disks 
have shown that they can generate sufficient en ergy to 
power a gamma-ray burst ('Popha m eTall I1999D . The 
viscosity in the disk pr imarily determines th e timescale 
of the resulting burst l)Naravan et al.l 12001(1 . while the 
mass of the BH may determine the overall energy scale 
(IPopham et alJ 119991) . These results are supported by 
ISetiawan et"an l|2004l I2005H. who find that the viscos- 
ity (and to a lesser extent, BH spin and disk mass) 
are responsible for determining the overall burst energy. 
Roughly speaking, luminosities of up to 10^° erg s~^ can 
be produced, so long as the disk is sufficiently massive 
(representing 5% of the total binary mass) and the ef- 
fective viscosity sufficiently strong {a = 0.1). 

Alternatively, general relativistic (GR) magnetohydro- 
dynamic (MHD) effects may allow infalling matter to tap 
the spin energy of a BH via the Blandford-Znajek effec t 
to produce an energetic jet l|Blandford fc Znaieklll977|) . 
In this case, the crucial parameters to determine the en- 
ergetics are the mass accretion rate and the spin of the 
BH. A revi ew of GRB mode ls tapping spin energy can 
be found in IZhang fc Meszar os (2004). 

2.1. Neutron Star-Neutron Star Binaries and 
Hypermassive Remnants 



There has been a great deal of numerical work per- 
formed to study merging NSNS binaries, including recent 
calculations performed using fully GR hydrodynamics 
and physically realis tic nuclear equations of state (EOS; 
iShibata et atll2005t iShibata fc Taniguchi) '20()6'l. From 
these results, two separate channels have been identified 
which could lead from a merger to the production of an 
SGRB. 

If the total binary mass is below some critical thresh- 
old, Mthr, a differentially rotating hypermassive neu- 
tron star (HMNS) can be formed that is initially stable 
against gravitational collapse, even though its mass is 
above the maximu m value f or uniformly rotating config- 
urations (Baumgar te et al] l2000: Morrison "6^01 12004). 
Delayed gravitational collapse follows, triggered either 
by gravitational wave dissipation if the remnant forms 
a bar ^ or by magnetic redistr i bution of angular momen- 



tum ( Baumgarte et al.l 12 000^ 'S hapirol 120001 iCook et alJ 
l200llDuez et al.ll2006al b: Shiba ta et al.ll2006(l . The nu- 



merical simulations of Shibata et al. use a stiff, realis- 
tic EOS (Akmal et al. 1998), for which the maximum 
mass is slightly higher than the highest meas ured pulsar 
mass M = 2.1 ±0.2 M© for PSR J0751-1807 l|Nice et alJ 
12005(1 . These simulations set a lower limit for the criti- 
cal mass Mthr > 2.7M0. This exceeds ah known NSNS 
binaries containing a radio pulsar, with the possible ex- 
ception of PSR 1913+16 (Mtot = 2.83 M©, Stairs 2004). 
Thus, it appears that HMNS formation is likely in most 
merging NSNS binaries. MHD simulations in full GR 
show that the HMNS undergoes a delayed collapse, re- 
sulting in a hot, magnetized torus surrounding a rotating 
BH, together with a magnetic field coUimated along the 
polar axis. These conditions are favorable for a burst 
powered either by ne utrino annihilation or MHD effects 
i(Duez et al.ll2006albl: IShibata et a'ni200^ . 

Alternatively, for higher-mass NSs, mergers with bi- 
nary mass ratios sufficiently far from unity might 
lead to the formation of a relatively massive disk 
around a BH formed pr omptly during the merger 
((Shibata fc Tanigucliil I2006() . For mass ratios q ~ 0.7, 
disk masses of 0.01 — O.IM© are possible, whereas more 
equal-mass mergers produce much lower mass disks, 
with insuffi cient thermal energy to power an SGRB 
((Shibata et al. 2005). 

The key unanswered question for the HMNS scenario is 
whether "baryon pollution" from merger ejecta along the 
polar axis can be cleared out of the funnel through which 
the presumed GRB j et will propaga te. Relativistic nu- 
meric al calculations l(Duez et al. 20 06albl : IShibata et alJ 
12006(1 indicate that MHD effects seem to be sufficient to 
accomplish this task, producing baryon densities along 
the polar axis t hat satisfy the constraints described in 
iShemi fc PiranI ((iggOl For the unequal-mass NSNS 
merger case, the polar axis is essentially free of interven- 
ing matter, but the likelihood of this scenario depends 
on how far from unity binary mass ratios are likely to 
fall, as all observed systems containing radio pulsars have 
q > 0.9. Systems with components too close together 
in mass produce significantly lower mass disks, which 
lack the required neutrino luminosity to power the SGRB 
(^hibata fc Taniguchi 2006). 



2.2. Black Hole-Neutron Star Binaries 
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The defining parameter for determining the quaUtative 
evolution of a BHNS merger is the binary mass ratio 
q = Mns/-^bh- Assuming that the tidal disruption of 
the NS begins at a separation or where its volume in 
isolation is equal to the volume of its Roche lobe, 

aR/MBH - 2.17^2/3(1 + qf /^C-' (1 ) 

for a NS of compactness C = M^a,/R^s l|Paczvnskill97H) . 
using geometrized units with G ^ c = 1. For a suffi- 
ciently large BH (and thus small q) is smaller than 
the innermost stable circular orbit (ISCO) aiscO: so that 
the NS passes through the ISCO before being disrupted 
tidally. For a typical NS of compactness C = 0.15 the 
critical mass ratio at which tidal disruption occurs at 
aisco is approximately q = 0.24. It has been argued that 
even for this "fiducial" binary the rapid plunge timescale 
at the ISCO might prevent significant disk formation 
l)Milleil l2005|) . Our simulations suggest, however, that 
rapid angular momentum transfer during tidal disruption 
ejects a significant fraction of the matter back outside 
the ISCO, leading to a sizable disk. Clearly, a detailed 
understanding of this dynamical process requires an ac- 
curate description of the strong field BHNS spacetime, 
as we discuss below. 

We performed 3-|-l-dimensional smoothed particle 
hydrodynamics (SPH) calculations of BHNS merg- 
ers in t he conform a l flat ness (CF) approximation 
to GR ijFaber et all l20M hereafter FBSTR). For 
a Schwarzschild BH, our scheme identifies including 
the relativistic ISCO exactly and accounts for rel- 
ativistic dynamics within the ISCO, unlike previous 
pseudo-New tonian calculations llLee & Kluzniak 1999; 
iJanka et all 119991: IRosswog et all 12004 ■Rosswoa,.2005^) . 
W e adopt the relativistic , irrotational binary models 
of iTaniguchi et alJ l)2005j) as initial data. Our adia- 
batic evolution scheme is the same as that described 
in FBSTR, but we solve the coupled non-linear ellip- 
tic field equations of the CF scheme using an iterative 
FFT-convolution based solver and add artificial viscos- 
ity to the hydrodynamical equations in order to study 
the shock heat ing of the material during the merger pro- 
cess (compare iHernaurst fc Katzill98^ . Both the initial 
data and the evolution scheme represent the state-of-the- 
art for including NS self-gravity self-consistently within 
a GR hydrodynamical formalism. 

We consider two models, both containing n — 1 poly- 
tropic NS. First is a NS of compactness 0.15, denoted 
"case A". As in FBSTR, our current method is lim- 
ited to extreme mass ratios, so we choose q = 0.1, for 
which Or = 3.2Mbh according to Eq. ^ We cannot di- 
rectly simulate the fiducial case {q — 0.24, C ~ 0.15) un- 
der these assumptions, so we take an alternate approach 
for our second calculation. Simultaneously reducing the 
mass ratio to g = 0.1 and the NS compactness to C = 0.09 
leaves flR/AfeH virtually unchanged. From the argument 
above, this binary, denoted "case B" should tidally dis- 
rupt slightly within the ISCO, mimicking the dynamics 
of the fiducial binary. 

The evolution of case A is straightforward and prob- 
ably uninteresting as an SGRB source: the entire NS 
spirals toward the BH, and no matter is ejected outside 
the ISCO to form an accretion disk. Case B, on the 
other hand, does lead to the formation of an accretion 
disk, as we show in Fig. ^ Note that radii are mea- 



sured in isotropic coordinates, in which the BH radius is 
r-BH = 0.5Mbh and the ISCO radius aisco ~ 5Mbh- We 
see that while the NS inspirals until 98% of its mass lies 
within the ISCO, rapid redistribution of angular momen- 
tum during tidal disruption causes an outwardly directed 
spiral arm to form, sending some matter back outside the 
ISCO. 

Eventually, the BH accretes Mace ~ 0.75AfNS directly, 
while part of the remaining mass forms a disk of mass 
Mdisk ~ 0.12Mns and part is ejected completely from the 
system {Mej ~ 0.13Mns)- Here we distinguish bound 
and unbound trajectories by the sign of uq — 1, where 
uo is the time-component of the matter 4- velocity, which 
remains nearly constant in time for outflowing gas on 
approximately geodesic trajectories. We note that this 
configuration satisfies all the geometric constraints re- 
quired for a GRB progenitor, as all matter lies in the 
equatorial plane rather than the polar axis. 

The bound matter in the disk is relatively cold at 
first, as the matter in the arm is initially ejected with- 
out strong shock heating. Over time, this disk generates 
heat via shocks as matter falls back and wraps around 
the BH forming a torus out to a radius of r 50AfBH 
within t = lOOOAfsH - .07 s; (see Fig. 12). The spe- 
cific internal energy in the inner part of the torus corre- 
sponds to a temperature T « 3 — lOAfeT^ w 2 — 7 x 
10^°if; (see Fig. I^J. The surface density in this re- 
gion is E « 2 — 3 X lO^^g c m~2. Assumi ng an opacity 
K = 7 X 10-i^(r/10"i4:)2 l|Di Matteo et al. 2002), we 
conclude that the disk is optically thick out to r ^ 15M. 
Using the diffusion limit, the neutrino flux is given by 
K. {7N^/3){(tT*/kT,) where = 3 is the num- 
ber of neutrino families, and a is the Stefan-Boltzmann 
constant. The neutrino luminosity is L^, w 2TTr'^F„ ^ 
2q54 gj.g g-i^ This value is roughly an order of magni- 
tude larger than that seen in the collapse of an HMNS 
l)Shibata et al. 2006), and should produce at least a com- 
parable annihilation luminosity, Lj^p ~ ^.O*^"^^ erg s~^. 

Qualitatively, the hot torus described here is similar 
both to the initial data used to study hydrodynamic disk 
evolu tion in earlier GRB models (Setiawan ct al. 20QJ, 
2005: iAlov et alJfeOOSj) . as well as to that for med by the 
collapse of an HMNS l|Shibata et alJ 120061: iDuez et alJ 
l2Q06b). However, we note that the torus described 
here is physically larger. Unlike previous pseudo- 
Newtonian calculations of BHNS mergers (J anka ct a^ 
IT999t iLee fc Khizniakl IT991 IRosswog et alJ 120041) . we 
find prompt disruption of the NS during the plunge 
(rather than a NS core which survives the initial mass 
loss phase) and a lower mass disk, but at a temperature 
similar to those calculati ons which include d a detailed 
microphysical treatment ijJanka et al.' 'l9991 . 

While we do not follow the long-term evolution of the 
accretion torus and surrounding material, we can esti- 
mate the fallback time for the bound component, as- 
suming geodesic orbits. Approximately 0.03A/ns should 
return back toward the BH on timescales equal to or 
longer than a second, which could in principle produce 
lower-energy bursts at later times. It is conceivable that 
this fallback accretion might explain the secondary X-ray 
flares observed in SGRBs many seconds after the initial 
burst (see Berger 2006 for a summary of the observa- 
tions), especially if self-gravity leads to the formation of 
higher density clumps of material, but further Simula- 



4 



J I I I I I I I I I I I I I I I M M I I I I I I I I I I I I I I I Mji n n — ' — I — ' — I — I — ' — 



- - 



— 5 - 




-5 K - 

iiiIiiiIiiiIiiiTIQ Ti I I I I I I I I I I I I I iT Zl I I I I I I L 



-20 - 



- 





02468 -4 -2 024 -20 20 

x/Mbh x/Mbh x/Mg^ 

Fig. 1. — Initial (left), early (center) and intermediate (right) configurations of the NS in case B, projected into the orbital plane. The 
NS has a compactness M/R = 0.09 and the binary mass ratio is Mns/A^bh = 0.1. We see the NS disrupts near the ISCO (dashed curve) 
to produce a single mass-transfer stream, which eventually wraps around the BH (solid curve) to form a torus. The initial orbital period 
is P = 105Mbh. 




Fig. 2. — Matter configuration at the end of the simulation, T = 990Mbh, projected onto the equatorial (left panel) and meridional 
(right panel) showing the hot torus located within r < SOMbh- Bound fiuid elements (satisfying uq — 1 < 0) are shown as crosses, unbound 
elements as points. Note the different scales. 



tions are required to establish this identification. 

Future observations should help to determine which 
progenitor candidates are responsible for the observed 
short GRB population. Gravitational wave measure- 
ments would provide important evidence if detected in 
coincidence: A GRB resulting from hypcrmassive col- 
lapse would occur noticeably delayed relative to the grav- 
itational wave signal from the inspiral and plunge phases 
of nearly equal-mass objects, whereas one resulting from 
a BHNS binary would occur almost immediately after 



the signal from a very unequal- mass merger. 
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Fig. 3. — Temperature T and surface density S as a function of cylindrical radius for the configuration shown in Fig. |^ The solid 
vertical line denotes the ISCO, and the dashed line the transition radius between matter optically thick to neutrinos within and optically 
thin outside. 
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